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We consider a system of two strongly coupled electron spins in zero magnetic field, each of which 
is interacting with an individual bath of nuclear spins via the hyperfine interaction. Applying the 
long spin approximation (LSA) introduced in Reff^ (here each bath is replaced by a single long 
spin), we numerically study the electron spin and entanglement dynamics. We demonstrate that 
the decoherence time is scaling with the bath size according to a power law. As expected, the 
decaying part of the dynamics decreases with increasing bath polarization. However, surprisingly it 
turns out that, under certain circumstances, combining quantum dots of different geometry to the 
double dot setup has a very similar effect on the magnitude of the spin decay. Finally, we show that 
even for a comparatively weak exchange coupling the electron spins can be fully entangled. 

PACS numbers: 76.20.+q, 76.60.Es, 85.35.Be 



I. INTRODUCTION 

The Loss-DiVincenco proposal is one of the most 
promising concepts for solid state quantum information 
processing. Here electron spins confined in semicon- 
ductor quantum dots are utilized as qubitsP^ The cen- 
tral drawback of this approach is the fast decoherence 
caused by the coupling of the electron spin qubits to 
the nuclea r spins o f the host material via the hypcrhne 
inter actior l 6 ' 7 ' 10 ' 1 5 -*^. For related reviews the reader is 
referred to RefsPSfSH. Other nanostructures in which 
similar situations arise are given by carbon nanotube 
quantum dots^, phosphorus donors in silicones! and ni- 
trogen vacancies in diamondP^HUl 

However, the hyperfine interaction allows to access the 
nuclear spins efficiently. Hence, when it comes to uti- 
lize them instead of the electron spins for quantum in- 
formation purposes, vice turns into virtue and the hy- 
perfine interaction gets a very advantageous character. 
Examples in this context are given by the possibil ity to 
built up an interface between light and nuclear spin d 30 ' 3 ^ , 
to polar ize nu clear spin baths^HII] to set up long-lived 
quantu m 1 35 ! 36 ' and classical memory devices or to gen- 
erate entanglement! 3 ^ 

Following the idea to take advantage of the hyperfine 
interaction, in a recent letter-^ we investigated a system 
of two exchange coupled electron spins, each of which 
is interacting with an individual bath of nuclear spins 
via the hyperfine interaction. In contrast to most of 
the appro aches c onsidered in the context of hyperfine 
interactio n 115 * 16 * 3 ^ ^, no magnetic field, enabling for a 
perturbative treatment of the problem, was applied to 
the electron spins. Using exact diagonalization studies, 
we demonstrated that the nuclear baths can be swapped 
and fully entangled, provided they are large enough. In 
order to be able to numerically consider the required 
system sizes, we introduced the so-called long spin ap- 
proximation (LSA). Here we assumed homogeneous cou- 
plings within each of the baths and considered them to 
be highly polarized. This allows to replace them by two 



single long spins. Interestingly, the spectrum of the two 
bath model with homogeneous couplings, studied in a 
preceding publication 4 ^, exhibits systematically degener- 
ate multiplets under certain conditions. Motivated by 
this, we distinguished between systems with and with- 
out inversion symmetry, i.e. a formal exchange of the 
central as well as the bath spins. In the latter case quan- 
tum dots of different geometry are combined to a double 
dot setup. Surprisingly, it turned out that here the swap 
performance is much better. 

In the present paper we apply the LSA in order to 
study the electron spin dynamics. The results comple- 
ment those of RefsP^H and in particular those of Ref. 43 , 
where we studied the electron spin evolution assuming 
the electrons to interact with a common bath of nuclear 
spins via homogeneous couplings. 

The paper is organized as follows: In Sec. |H]we intro- 
duce the model and the methods. In particular, we in 
detail discuss the applicability of the LSA with respect 
to the electron spin dynamics. We then study the spin 
and entanglement dynamics in the limit of an exchange 
coupling which is much larger than the hyperfine energy 
scale. Here the nuclear baths act as a perturbation. This 
is a particularly interesting case, as exceptionally long 
decoherence times can be expected. In Sec. |III| we fo- 
cus on the time evolution of the electron spins. In a 
first step we study basic dynamical properties. In par- 
ticular we demonstrate that in certain parameter ranges 
the process of decoherence is incomplete. Furthermore, 
we find a simple empirical rule describing the dynamical 
signatures of different initial states. We then quantita- 
tively investigate the decoherence time and the magni- 
tude of the spin decay. As expected from Ref.^21, the 
decoherence time scales with the system size according 
to a power law. As already mentioned, in Ref!^ it was 
demonstrated that the nuclear spin dynamics strongly 
benefits from combining quantum dots of different geom- 
etry to the double dot setup. In full generality, this result 
can be confirmed only in certain parameter regimes. In 
Sec. |IV| we then focus on the entanglement dynamics and 
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FIG. 1: Illustration of the one-bath and the two-bath model 



demonstrate that, surprisingly, even for a comparatively 
weak exchange coupling the electron spins can be fully 
entangled. 



II. MODEL AND METHODS 

The hyperfme interaction in a double quantum dot is 
described by the Hamiltonian 



H 
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Ii + JexSl ■ S2, (1) 



where S!j are the electron and Ii are the nuclear spins. 
The parameter J ex denotes an exchange coupling be- 
tween the two electron spins, which can be adjusted in 
a range of [-1CT 3 , 10~ 3 ]eV. The constants Af ] are 
the hyperfme couplings of the two electron spins. In a 
realistic quantum dot, these are proportional to the elec- 
tronic wave function of the j-th electron at the site f*j of 
the i-th nuclear spin: 
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For typical GaAs quantum dots this leads to an interac- 
tion with N ~ 10 6 nuclear spins and the overall hyperfme 
coupling strength of the j-th electron, 
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is of the order of [lCT 4 , l(T 5 ]eV (see RefP). 

Due to the spatial variation of the electronic wave func- 
tion, the hyperfme couplings are clearly spatially depen- 
dent. However, for any set of hyperfme coupling con- 
stants the Hamiltonian, obviously, conserves the total 
spin J = Si + S2 + X)i=i Ii- This is a very helpful sym- 
metry for exa ct nu merical diagonalizations of the Hamil- 
tonian matrbPSHiJ through which we will gain the dy- 
namics of the system in what follows. Here we consider 
the eigensystem of the Hamiltonian 
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and decompose the initial state \a) into a sum of energy 
eigenstates: 



(5) 



Applying the time evolution operator U — e~R and 
tracing out the nuclear degrees of freedom then gives the 
reduced density matrix for the electrons, 
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from which the dynamics of all observables can be calcu- 
lated. For further details see RefP2J 

There we investigated the case of two electron spins 
coupled to a common nuclear spin bath. In what follows, 
however, we consider the case of two separate baths as 
depicted schematically in Fig. [l] In the first case the 
two electron spins are assumed to be very close to each 
other so that both interact with the same group of nuclear 
spins, whereas in the present case they are spatially more 
separated, leading to an interaction with an individual 
group. The realistic situation of a double quantum dot 
will of course lie between these two extreme cases. 

As already mentioned, in the present paper we apply 
the LSA to the two bath system. In the following sub- 
sections we give a detailed discussion of the model with 
a particular focus on its limitations. 



A. The long spin approximation (LSA) 

Let us consider two separate spin baths of equal size 
with homogeneous couplings to one of the two electron 
spins each and introduce Ij = J^i=i ^ij 1 where the 7y are 
the Nj nuclear spins the j-th electron spin interacts with. 

This means that N = ATj + N 2 and J2iLi h = h+h, 
where, for simplicity, we will consider N\ = N 2 in what 
follows. Now the squares of the total spin of each bath 
are separate conserved quantities. Moreover, the same 
holds for the square of any sum over a subset of spins of 
each bath, 
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where we have, for the sake of brevity, denoted the set 
of all the latter operators of the j-th bath as {Kj }. The 
corresponding quantum numbers {Kj} can be used to 
characterize specific Clebsch-Gordan decompositions of 
each bath. 

The initial state \a) is given by a direct product be- 
tween the initial state of the electron spins \a e ) and the 
initial state of the baths \a n ). Provided the two dots 
are spatially well-separated, the two resulting baths have 
to be considered as practically uncorrelated. Hence, the 
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caused by the flip-flop terms 
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in H or -S/lsa, respectively. This is exactly the part of 
the Hamiltonian, which is eliminated in most of the ap- 
proaches by applying a s trong magn etic field to t he cen- 
tral spin system (see RefsPEMUS]) In Ref s EIlH 

we also 

concentrated on the dynamics which are purely due to the 
flip-flop terms. 



FIG. 2: The two baths are approximated by two long spins. 



B. Homogeneous couplings on long times scales 



state of the nuclear baths is again a direct product be- 
tween the states of the two baths, \/3j). In general such 
a state reads 
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where \Ij,rrij, {Kj}) are the eigenstates of Ij. If the re- 
spective bath is now strongly polarized, the number of 
contributing multiplets in Q drastically decreases.^ If 
we are close to full positive or negative polarization, we 
can drop the quantum numbers {Kj} and consider the 
initial state to be given by \I, rrij) with rrij ~ ±1. Due 
to Q , no "cross terms" between different multiplets con- 
tribute to the dynamics and all physics is then captured 
in the LSA Hamiltonian 
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S 2 -I 2 + J ex SfS 2 , (9) 



sketched in Fig. [2] The coupling constants A^ /2I re- 
sult directly from ^ by considering 7y = 1/2: As all 

couplings A^p are chosen to be equal to each other, |s| 

yields A^ = A^ /Nj. The quantum number Ij ranges 
from to Nj/2. As our model is based on highly polar- 
ized baths, we choose the maximal value. Together with 
Ni = N 2 this yields A\ = A® /2I. Although the LSA 
Hamiltonian is not exactly solvable, the approximation of 
the baths by single long spins reduces the dimension of 
the problem so that exact numerical diagonalizations are 
possible on arbitrary subspaces even for comparatively 
large baths. 

Accounting for the J z symmetry, the nuclear state \a n ) 
explicitly reads 



\a n ) = \I,M -m-i)\I,i). 



(10) 



Here M denotes the quantum number associated with 
J z and m the one related to S z . The parameter i is 
introduced in order to account for the deviations from 
rrij = ±1. Hence, it has to be chosen in the vicinity of I 
or (M — m — I), respectively. Note that for an initial state 
which is a simple product state like (10), all dynamics is 



In Ref ESI we considered the one-bath model illustrated 
in the upper panel of Fig. [T] for homogeneous couplings 
and initial states with a very low bath polarization := 
(N - 2N D )/N of (l/N). Here N D denotes the number 
of flipped spins in the bath. The central spin dynamics 
shows periodic behaviour. This clearly has to be regarded 
as an artifact caused by the homogeneity of the couplings. 
For short time scales, meaning times much smaller than 
the recurrence time, the results for decoherence times 
found there compare well with experimental values. 

As explained above, within the LSA we assume the 
couplings to be homogeneous and the baths to be highly 
polarized. However, high polarizations naturally lead to 
long time scales for the electron spin decoherence times. 
Consequently, it has to be analyzed to what extent the 
two assumptions of the LSA contradict each other. In 
Ref.El we already investigated this question with respect 
to the nuclear spin dynamics, where we considered a 
Gaudin model, as corresponding to one of the first terms 
in ([!]). We found that, qualitatively, inhomogeneities be- 
come less important with increasing polarization. Typi- 
cally, in such a context one would give a quantitative ar- 
gument by evaluating the fidelity (to be precisely defined 
below) rather than studying the dynamics on a quali- 
tative level. However, with respect to the nuclear spin 
dynamics considered in RefP this does not make sense, 
obviously, as the bath consists of many spins so that a 
certain value of (I z ) can be realized by a whole set of 
nuclear states. 

In the following, we again consider a usual Gaudin 
model and investigate the time-averaged fidelity F with 
respect to homogeneous and inhomogeneous couplings 
via exact diagonalization. This is given by 



F =hl dt\(U h a\U ih a)\ 

1 Jt=Q 



(12) 



with Uh and C/jh being the time evolution operators for 
the homogeneous or inhomogeneous Hamiltonian respec- 
tively. We choose an initial state which is a direct product 
between an electron spin pointing upwards and a ran- 
domly correlated bath state. This is a superposition of 
all possible states with (in our case) real random coeffi- 
cients, which we choose in the interval [—1, 1]. Randomly 
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correlated states lead to highly repr oduci ble results and 
can therefore be regarded as generic! 20 * 44 ^ 

The results are shown in Fig. [3] We fix three different 
system sizes of N — 12,20,30 bath spins for a reason- 
ably long period T — 400(h/A). In the left panel we 
plot F against In order to get a better comparison 
between the different system sizes, in the right panel we 
show the same data plotted against the bath polarization 
Pb- Obviously, the fidelity is strongly increasing with the 
bath polarization. Furthermore, it decreases with an in- 
creasing number of bath spins. Note that the different 
curves approach each other with increasing number of 
bath spins. 

The highest experimentally feasible polarizations are 
around 80%, as reported in RefPM. On first sight, the 
results shown in Fig. [3] indicate that even for such 
high polarizations considering homogeneous couplings on 
comparatively long time scales is restricted to extremely 
small systems. This would strongly contradict the pur- 
pose of the LSA. However, it turns out that the fidelity 
is an extremely sensitive measure underestimating the 
applicability of the LSA: In Fig. [4] we plot the spin dy- 
namics for inhomogeneous and homogeneous couplings 
for Nd = 4 pb = (1/3), corresponding to a very low 
fidelity of F = 0.055194. Such a small value clearly sug- 
gests that the dynamics in the inhomogeneous and the 
homogeneous case are fundamentally different. The am- 
plitude of (S z (t)) decaying to zero without any recurrence 
on the considered time-scales would be an example of a 
natural expectation for the first case (compare e. g. th e 
results presented in RefsPHU with those of RefePESJ). 
Furthermore, one would guess that the time-averaged 
values of (S z (t)} in the inhomogeneous and the homo- 
geneous case strongly differ from each other. However, 
as can been seen from Fig. 4, neither of these expecta- 
tions are met. This means that even very small fidelities 
correspond to a rather good qualitative agreement of the 
dynamics. Considering highly polarized baths, as done 
within the LSA, it is therefore justified to choose homoge- 
neous couplings even on comparatively long time scales. 



III. ELECTRON SPIN DYNAMICS 

In the following we restrict ourselves to the limit 
J ex /(A/2I) > 1, where we defined A = A^+A (2 \ Here 
the dynamics is dominated by the electron spin coupling 
term and the baths act as a perturbation. As a conse- 
quence, long decoherence times, enabling e.g. to fully 
entangle the two electron spins, have to be expected. 
We will distinguish between a "strong coupling" and an 
"ultra strong coupling" limit. In the first case an only 
moderately large exchange coupling, J ex /A « 1, is con- 
sidered so that the condition J ex /(A/2I) ^> 1 is realized 
mainly through the length of the bath spins, whereas in 
the second case we choose a very strong exchange cou- 
pling, meaning that here we already have J ex /A 3> 1. As 
It = I2, a zero "detuning" A :— A 2 — A 1 is associated 
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FIG. 3: The time-averaged fidelity is plotted against the 
number of flipped spins in the bath (left panel) and against 
the bath magnetization (right panel). We consider a ran- 
domly correlated initial bath state with coefficients in [—1, 1] 
for TV = 12, 20, 30. The fidelity strongly increases with in- 
creasing polarization. Furthermore, it decreases with an in- 
creasing number of bath spins, where the decrease gets the 
weaker the larger the baths become. 
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FIG. 4: Central spin dynamics for inhomogeneous and ho- 
mogeneous couplings for TV = 12 and TVd = 4 <S> pt = (1/3), 
as corresponding to F = 0.055194. Just as in Fig. [3] we con- 
sider a randomly correlated initial bath state with coefficients 
in [—1,1]. In both cases the spin is oscillating around a very 
similar mean value, where the amplitude is significantly larger 
for the case of homogeneous couplings. Hence, even in this 
extreme case of a very small fidelity, the dynamics are still 
somewhat similar to each other. 
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FIG. 5: Electron spin dynamics for the strong and the ultra 
strong coupling case J ex / (A/21), J BX / A — 350 with A = 
in the left column and A 7^ in the right one. We consider 
I = 120 and M r = 0.17. The initial state is given by \a) = 
)\M — I, I). As illustrated by the dotted lines, we choose the 
point at which the amplitude does not change anymore as the 
decoherence time. We find very regular dynamics, where the 
decoherence times are obviously larger in the case of broken 
inversion symmetry A / 0. Furthermore, the oscillations do 
not fully decay in the ultra strong coupling case. 



with a system invariant under inversions 1 o 2 P Phys- 
ically a detuning different from zero corresponds to dots 
of different geometry combined to a double quantum dot. 
Throughout the paper we will consider the dynamics on 
subspaces with fixed magnetization. It therefore suffices 
to investigate the z components of the spins. Further- 
more, due to their strong coupling, the dynamics of the 
two electron spins can be read off from each other even 
for A ^ 0. Therefore, we always focus on the time evo- 
lution of the first electron spin. Note that the energy 
scale is given by (A/21) and consequently the time will 
be given in units of hi /it A. 



A. Basic dynamical properties 

In order to give a basic impression of the dynam- 
ics, in Figs. [5] and |6] we fix J ex /(A/2I) = 350 for the 
strong and 3 ex /A = 350 for the ultra strong coupling 
case and plot the dynamics of the first electron spin for 
/ = 120. We consider the relatively low "magnetization" 
M r := M/(2I + 1) = 0.17. This means that we concen- 
trate on initially nearly antiparallel baths. We study the 
inversion invariant case A = as well as A ^ 0. All ini- 
tial states considered in Figs. [5]and[6]have an antiparallel 
electron spin configuration \a e ) — In Fig. [5]we con- 

sider an initial nuclear state with a maximally negative 
z component of the first bath spin, Ii . This corresponds 



to i = I in ( 10 1. We clearly see that in the ultra strong 



FIG. 6: Electron spin dynamics for the strong coupling case 
J ex /(A/2I) = 350 and A = 0.015. We consider I = 120 
and M r — 0.17. The initial states are given by \a) = |"f|\|J. 
)\M - I + i, I - i) with i = 1,2, 3, 4. We find dynamics with 
an envelope, decaying in a quite similar way to the one for 
i = shown in Fig. [5] However, here additional beatings 
occur. Their number is equal to i. 



coupling limit the time evolution for initial states of the 
above mentioned form does not fully decay. Indeed, for 
small, inversion invariant systems this is also the case in 
the strong coupling limit. Varying i slightly away from J, 
the dynamics in the ultra strong coupling case does not 
show any qualitative change. As can be seen in Fig. |6j 
this is also the case for the envelope of the dynamics in 
the strong coupling limit. However, here additional beat- 
ings occur. Surprisingly, there is a clear empirical rule 
concerning these additional low frequency oscillations: If 
the z component of the first bath spin deviates by i from 
the maximal negative value, \a n ) = \M — I + i, I — i), 
the dynamics shows exactly i beatings. In Fig. [6] the 
case of broken inversion symmetry is considered, where 
the beatings are particularly pronounced. At the time 
being, we are not able to explain this effect. However, 
it seems that non-trivial dynamical regularities are typi- 
cal for central spin models with homogeneous couplings. 
Indeed, in Refp2l we reported on a rule for the one bath 
model, which relates the number of flipped spins in the 
initial state of the bath to the number of local extrema 
in the oscillations of the central spins. Also the dynamics 
has been calculated on a fully analytical level, we have 
not been able to give an explanation of these regularities. 



B. Decoherence time and magnitude of the spin 
decay 

In direct analogy to the investigations in RefP^, in the 
following we investigate the scaling of the decoherence 
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FIG. 7: Scaling of the decoherence time with the bath spin 
length / on a double logarithmic scale for the strong cou- 
pling case. Different coupling ratios J ex / {A/21) and detun- 
ings A are considered. We choose the initial state |a) = |"f|\|J. 
)\M - 1,1) for M r = 0.17. We find power laws ~ with 
v = 1.29 (Jen/ {A/21) = 350, A = 0), v = 1.03 {J ex /{A/2I) = 
175, A = 0), v = 1.54 {J ex /{A/2I) = 350, A = 0.015A), and 
v = 1.65 {J ex /{A/2I) = 350, A = 0.019A). Although the 
limit Jex/{A/2I) 3> 1 is left unaltered, surprisingly the scal- 
ing changes with the coupling ratio. Breaking the inversion 
symmetry leads to an increase of the parameter v. 



FIG. 8: Scaling of the decoherence time with the bath spin 
length J on a double logarithmic scale for the ultra strong 
coupling case. We fix J ex /A = 350 and vary A weakly away 
from zero. We choose the initial state \a) = \ili}-)\M — I , I) for 
M r — 0.17. The curves are again fitted to power laws ~ I~ u 
with v = 0.79 (A = 0), v = 0.49 (A = 0.00L4), and v = 0.54 
(A = 0.006^4). Note that the absolute of the decoherence 
times are smaller than in the strong coupling case because 
the dynamics does not fully decay. Breaking the inversion 
symmetry leads to a decrease of v to v w 0.5, which is the 
value found in Ref 43 for the one bath model. 



time with the spin length. It is clear that such an in- 
vestigation can not yield perfectly reliable values, as the 
spin length is of course restricted to comparatively small 
values due to the limited computational power. Con- 
sider for example M r — 0. Here the dimension of the 
Hilbert space is given by (8/ + 2), limiting the length of 
the spins to values of the order of I ~ 10 2 . Still, the re- 
sults give a clear idea about the type of scaling and allow 
for a qualitative comparison between different parameter 
regimes. In the following we concentrate on initial states 
| a) = IM) \M-I, I) for M r = 0.17. As already explained, 
i has to be in the vicinity of I and the envelope remains 
unaffected when varying i slightly away from its maximal 
value. Hence, the results for \a) = \^ij.)\M — I , I) can be 
regarded as generic. As can be seen from Fig. [5j in the 
ultra strong coupling limit (Sf(t)} does not decay to a 
constant value, but oscillations of quite regular shape 
remain, i.e. the decoherence process is not complete. 
Therefore we define the time from which on the ampli- 
tude does not change anymore as the decoherence time. 
Numerically this is realized by dividing the time axis in 
intervals with a length larger than the period of the reg- 
ular oscillations and determining the maximal value in 
each interval. If this value does not change anymore over 
a fixed number of intervals, the lower bound of the first 
interval in which the respective value appeared is chosen 
as the decoherence time. In the left panels of Fig. [5] this 
choice is illustrated by the dotted lines. 



In Ref. 43 it has been shown for the one bath model 
that the decoherence time scales with the size of the 
bath according to a power law ~ N~ v . Indeed, we 
find the same behavior for the present case. In Figs. 
[7J [8] the decoherence times for the strong and the ultra 
strong coupling case are plotted against the spin length 
J on a double logarithmic scale. For the inversion sym- 
metric case we consider J ex /(A/2I),J ex /A = 350 and 
J ex /(A/2I) = 175, where it obviously does not make any 
sense to choose a second value for the ultra strong cou- 
pling limit. For the case of broken inversion symmetry we 
fix J ex /(A/2I), J ex /A = 350 and fix A = 0.015A, 0.019A 
for the strong coupling and A = 0.001A, 0.006^4 for the 
ultra strong coupling case. The values for the latter are 
chosen to be particularly small, because, as exemplified 
in the bottom panels of Fig. [5] the dynamics is highly 
sensitive with respect to a change of the detuning and 
becomes completely coherent on any relevant time scale 
for larger values. 

As expected, for the strong coupling case the deco- 
herence time is scaling much stronger than for the ultra 
strong coupling limit. Note that the values for the ultra 
strong case are much smaller only due to the fact that 
here the dynamics does not fully decay. As can be seen 
from Fig. [7| in the strong coupling limit the scaling does 
change significantly with the coupling ratio J ex /(A/2I). 
This is surprising as a small change in the ratio leaves the 
limit J ex /(A/2I) 3> 1 unaltered and hence one would ex- 
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FIG. 9: Scaling of the amplitude, given by the difference 
of a local maximum and the following local minimum, with 
the magnetization and the detuning for the cases, where 
the dynamics does not fully decay. In the left panel we fix 
A = 0. For the strong coupling case we consider I — 20 and 
J ex /(A/2I) = 350, 175. For the ultra strong coupling case 
we fix I = 80 and J ex /A — 350. We choose the initial state 
a ) ~ 1^4)1-^" — I, I)- In all cases we find a linear dependence 
of the amplitude on the magnetization with gradients a close 
to one. The values are given a = 0.97 (J ex /(A/2I) = 350), 
a = 0.92 (J ex /(A/2I) = 175), and a = 1.0 (J ex /A = 350). In 
the right panel we plot the the amplitude against the detuning 
for M r w 0.17. Here we find a highly non-linear dependence. 



responds to nuclear spin baths highly polarized in either 
the same or opposite directions. However, in Fig. [9] we 
plot the amplitude, defined as the difference between a 
local maximum and the following local minimum, for the 
whole range of M r , which corresponds to either paral- 
lel or anti. We consider the strong as well as the ultra 
strong coupling case. For the first one we fix I = 20 and 
two values for the coupling ratio J ex /(A/2I) = 350, 175. 
In both cases we find a linear dependence with a gra- 
dient close to one, meaning that the ratio does not sig- 
nificantly influence the decaying part. Concerning the 
ultra strong coupling case, we set I = 80 and consider 
Jex/A = 350. The scaling is practically identical to the 
one for the strong coupling case. As already discussed in 
the preceding section, we found that for the ultra strong 
coupling case the decaying part is not only influenced by 
the magnetization but also by the detuning. In the right 
panel of Fig. [9] we plot the amplitude against the detun- 
ing for a fixed magnetization M r — 0.17. In contrast to 
a variation of the magnetization, here we find a highly 
non-linear dependence, which can not be fitted by some 
simple power law. 

In RefP we demonstrated that a non-zero detuning is 
very advantageous with respect to swapping and entan- 
gling the nuclear spin baths. When it comes to the elec- 
tron spin dynamics, however, in general this is the case 
only the ultra strong coupling limit. 

IV. ENTANGLEMENT DYNAMICS 



pect the scaling to be insensitive against a change of the 
coupling ratio. Furthermore, the absolute values of the 
decoherence time clearly decrease with decreasing cou- 
pling ratio J ex /{A/2I) as expected. However, as a coun- 
terintuitive effect the scaling with the system size turns 
out to be weaker for the smaller of the two ratios. Break- 
ing the inversion symmetry has a significant effect in the 
strong as well as the ultra strong coupling limit. In the 
first case the exponent v increases, whereas in the latter 
it decreases to v 0.5. This is the value derived in RefJ™ 
for the one bath model. 

As explained in the preceding subsection, in the ul- 
tra strong coupling case the dynamics does not show full 
decoherence. If the spin length / is small and we have 
A = this is also the case for only strongly coupled elec- 
tron spins. We now analyze the scaling of the decaying 
part of the dynamics as a function of the magnetization 
and, in the ultra strong coupling case, as a function of the 
detuning A. Concerning the strong coupling case the re- 
sults are, obviously, only of fundamental interest. Even 
in SiGe and carbon based quantum dots (only around 
4.7% of the Si isotopes are spin carrying™, in carbon 
even only 1%E^), the electron spins interact with a few 
thousands of nuclear spins. 

Note that our LSA model is valid only for relatively 
small and relatively large magnetizations M r , which cor- 



We now close the discussion of the electron spin dy- 
namics with an investigation of the entanglement be- 
tween the two electron spins. In order to quantify the 
non-classical correlations, we consider the concurrence 
defined bjESl 



C(t) = max 




where Xi are the eigenvalues of the non-hermitian ma- 
trix p e (t)p e (t) in decreasing order. Here p e (t) is given by 
(tTj, ® a y ) p* e {t) (<J y ® a y ), where p* e {t) denotes the com- 
plex conjugate of p e (t)-the reduced density matrix of the 
electrons as defined in 

In the following, we ask to what extent it is possible 
to entangle initially uncorrelated electron spins. There- 
fore, we again consider initial states with electron spin 
configurations |flJJ-). In particular we are interested in 
a lower bound for the ratio J ex /(A/2I), meaning that 
we adjust the couplings J ex , A±, A2 to the lowest possi- 
ble ratio so that the concurrence still becomes equal to 
one. As to be expected, the lower bound lies in the ultra 
strong coupling limit. However, surprisingly it turns out 
that it is not determined by the ratio J ex /A, but only by 
J ex /max{Ai, A-{\. The concrete value of this ratio de- 
pends on the initial state of the nuclear spins. An upper 
bound is given by the (,as explained above, unphysical) 
case of randomly correlated states. As an empirical rule 
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size of the baths according to a power law. As expected, 
it turned out that the decaying part decreases with in- 
creasing polarization. However, surprisingly it also de- 
creases with increasing detuning, provided the electrons 
are bound ultra strongly. Hence, with respect to the 
electron spin dynamics the advantageous character of a 
non-zero detuning, found in RefP^ for the time evolution 
of the nuclear baths, can only be confirmed in the ultra 
strong coupling limit. Finally, we demonstrated that it 
is possible to fully entangle the electron spins even for a 
comparatively weak exchange coupling. 
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FIG. 10: Time evolution of the concurrence between the 
two central spins for I = 40. We consider the case of a 
low polarization (M r — 0.09) in the left panel and the case 
of a high polarization (M r = 0.86) in the right panel. In 
both cases we choose a detuning of A = 0.019^1 and consider 
Jex/A 2 = 35/4, which satisfies {I!} , and J ex /A 2 = 10/4, 
which violates (141. The nuclear spins are in a randomly cor- 
related state with coefficients in [—1, 1] initially. One clearly 
sees that if (141 is satisfied, the concurrence becomes one, 
whereas for a stronger coupling to the baths the electron spins 
can not be fully entangled. 



of thumb here we find: 



max{^!, A 2 } 



> 8.8 (14) 



In Fig. 10 we illustrate the rule by plotting the dynamics 
for randomly correlated initial states with coefficients in 
[—1,1] by considering parameters satisfying and violat- 



ing ( 14 ) . We choose a rather small system of / = 40 and 
concentrate on the case of broken inversion symmetry 
A = 0.019A We plot the time evolution for a low po- 
larization of M r — 0.09 in the left panel and fix a rather 
high polarization of M r — 0.86 in the right panel. It is 
visible that the maximal value of the entanglement drops 



slightly under one if ( 14 ) is violated 



V. CONCLUSION 



In summary, we numerically studied the electron spin 
and entanglement dynamics in a system of two strongly 
coupled electron spins, each of which is interacting with 
an individual bath of nuclear spins via the hyperfine in- 
teraction. We applied the LSA introduced in RefPJ where 
the two baths are replaced by two single long spins, and 
focused on the limit of an exchange coupling much larger 
than the hyperfine energy scale. Here we distinguished 
between a strong and an ultra strong coupling case. We 
demonstrated that the decoherence time scales with the 
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